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modifications  to  account  for  the  frequency  dependent 
properties  of  viscoelastic  materials  (loss  factors  and  shear 
modulus)  using  MSC/NASTRAN.  A  modal  strain  energy 
approximation  method  was  also  applied  to  each  plate. 
Comparisons  of  the  frequency  response  analysis  to  the  modal 
strain  energy  method  were  made  for  each  plate  structure. 
The  experiments  were  performed  on  the  plates  to  establish 
the  base  line  data  for  the  comparisons. 

The  modal  strain  energy  method  frequency  response  plots 
showed  good  agreement  with  the  frequency  response  analysis 
method  plots  for  all  plate  structures  analyzed.  Plate 
structures  with  fewer  natural  frequencies  in  the  frequency 
range  investigated  showed  better  agreement  between  the  two 
prediction  methods  than  the  structures  with  more  natural 
frequencies.  The  experimental  data  for  the  plate  analyzed 
indicate  less  damping  than  calculated  by  either  of  the 
finite  element  prediction  methods.  The  experimentally 
determined  natural  frequency  modes  were  well  approximated  by 
both  of  the  finite  element  prediction  methods. 
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I.   INTRODUCTION 

Advances  in  the  material  sciences  area  have  resulted  in 
metals  with  high  strength  to  weight  ratios.  The  application 
of  such  material  can  result  in  structures  with  low  levels  of 
damping.  Exposure  of  these  structures  to  high  levels  of 
excitation  creates  a  need  for  damping  techniques  which  are 
both  inexpensive  and  effective.  Introduction  of  a 
distributed  viscoelastic  material  with  an  incorporated 
constraining  layer  has  gained  widespread  acceptance  as  a 
means  of  reducing  vibration  response.  Many  finite  element 
codes  have  been  developed  which  incorporate  frequency 
prediction  methods  for  use  in  design  applications,  reducing 
the  need  for  experimental  analysis.  However,  accurate 
behavioral  predictions  of  applied  viscoelastic  material  have 
remained  costly,  due  to  the  frequency  and  temperature 
dependence  in  viscoelastic  materials. 

The  response  of  the  simplest  structures  can  be 
represented  by  continuous  classical  methods  but  the  finite 
element  approach  remains  the  method  of  choice  when  complex 
physical  systems  must  be  analyzed.  The  use  of  a  finite 
element  method  to  find  the  system  response  requires 
modifications  to  capture  the  environmental  dependencies  of 
the  applied  viscoelastic  layer.  One  standard  method  of 
conducting  a  frequency  analysis  of  dampened  structures  is 


through  repeated  application  of  the  frequency  response 
analysis  over  incremented  steps  in  frequency.  A  second 
approach,  using  modal  strain  energy  approximations  in 
conjunction  with  the  finite  element  method,  has  proved 
effective  in  predicting  excitation  responses  of  structures 
with  a  constrained  viscoelastic  layer.  This  paper 
investigates  the  two  prediction  methods  described  above  as 
applied  to  plate  structures.  MSC/NASTRAN  was  chosen  for  the 
analyses  because  of  its  widespread  availability  and  previous 
use  in  this  area.  The  prediction  methods  were  applied  to 
two  plate  structures,  one  with  the  damping  material  having 
edge  exposure  to  the  environment  and  a  similar  plate  with 
the  damping  material  protected  from  the  environment  and 
possible  corrosion,  through  the  use  of  an  edge  weld.  The 
edge-welded  plate  was  evaluated  before  and  after  welding  to 
determine  possible  damage  to  the  viscoelastic  layer.  The 
prediction  methods  were  evaluated  for  ease  of  application 
and  prediction  accuracy  as  compared  to  the  experimental 
data . 


II.   THEORY  AND  BACKGROUND 

Continuous  structures  in  the  physical  world  possess 
distributed  characteristics  of  mass,  damping  and  stiffness. 
[Ref.  l:pp.  140-141]  Structures  with  such  characteristics 
are  said  to  contain  an  infinite  number  of  degrees  of 
freedom.  The  classical  vibration  analysis  of  such  a 
continuous  system  involves  the  construction  of  a  discrete 
mathematical  model  with  a  finite  number  of  parameters  to 
approximate  the  overall  system  response.  Construction  of  a 
model  of  N  discrete  elements  allows  the  discretized 
continuous  system  to  be  represented  by  the  coupled  matrix 
equations: 

CM]{q(t)}  +  [C]{q(t))  +  [K]{q(t)}  =  {Q(t)}         (2.1) 

Linear  transformation  and  normalization  of  the  undamped 
portion  of  Equation  (2.1)  and  assumption  of  damping 
approximated  by  an  appropriate  diagonal  matrix  gives  the 
uncoupled  equations: 


(a(t)}  +  [  2^ajn  ]  (  o'(t)  )  +  [  oo  n  ]  (  a(t)  }  =  {l(t))      (2.2) 


{q(t)}  =  [(})]{  a(t)}  (2.3) 


where : 

{q(t) }  =  displacement  vector 

[  (j)  ]     =  mode  shape  matrix 

[  (1)^  ]  =  vector  of  natural  frequencies 
n 

C      =  damping  factor 

{l(t)}  =  transformed  matrix  of  applied  loads. 

Application  of  modal  analysis  through  developed  finite 
element  codes  have  resulted  in  readily  obtainable  solutions 
to  Equations  (2.2)  and  (2.3)  for  complicated,  lightly  damped 
structures.  [Ref.  2: Sees.  2.3,4.3.2;  Ref.  3:p.  2] 
Recently,  viscoelastic  materials  have  been  added  to 
structures  to  improve  their  overall  damping  characteristics. 
However,  the  addition  of  these  viscoelastic  materials  to 
design  structures  introduce  both  frequency  and  temperature 
dependent  variables  that  must  be  incorporated  into  the 
finite  element  analysis.   [Ref.  4:p.  2] 

Viscoelastic  damping  has  been  observed  in  many  polymeric 
materials  and  has  been  directly  attributed  to  the  relaxation 
and  recovery  of  long  carbon  atom  chains  of  the  polymer  after 
deforming.  [Ref.  5:p.  67]  Because  of  their  effect  on 
molecular  motion,  temperature  and  frequency  influence  these 
"relaxation  and  recovery"  properties.  Thus,  a  complete 
understanding  of  the  temperature  and  frequency  dependent 
behavior  of  the  material  is  essential. 


Moderate  changes  in  temperature  or  large  changes  in 
frequency  result  in  definite  predictable  alterations  in  the 
viscoelastic  material  properties  of  shear  modulus  of 
elasticity  and  loss  factor.  Actual  quantitative  changes  in 
the  damping  properties  vary  substantially  with  the  chosen 
material.  Qualitative  changes  are  categorized  into  three 
descriptive  regions;  the  glassy  regions,  transition  region, 
and  the  rubbery  region  as  shown  in  Figure  2.1.  In  a  given 
viscoelastic  material,  all  three  regions  are  typically 
contained  in  two  to  three  decades  of  frequency.  [Ref.  5: p. 
91] 

Damping  at  relatively  low  temperatures  for  the  given 
material  is  in  the  glassy  region  and  exhibits  a  low  loss 
factor  and  high  shear  modulus  of  elasticity.  As  the 
material  temperature  is  increased  the  shear  modulus  begins 
to  drop  slowly  as  the  loss  factor  rises.  Sufficient 
temperature  increases  allow  the  damping  material  to  enter 
the  transition  region,  characterized  by  a  rapidly  decreasing 
shear  modulus  value.  The  loss  factor  reaches  a  peak  value 
at  the  maximum  rate  of  change  in  shear  modulus  value. 
Further  increases  in  temperature  result  in  the  material 
entering  the  rubbery  region  with  a  corresponding  decrease  in 
loss  factor  and  shear  modulus  values.  While  the  above 
discussion  is  descriptive  of  varying  properties  with 
temperature  increases,  it  also  applies  to  frequency 
decreases.    Examination  of  several  materials  shows  a  25 


degree  fahrenheit  increase  being  roughly  equivalent  to  a  one 
decade  frequency  decrease.  For  example,  a  25°  F  increased 
temperature  is  roughly  equivalent  to  a  10  kHz  to  1  kHz  or  a 
1000  Hz-lOO  Hz  drop  in  frequency.   [Ref.  6:pp.  A1-A21] 

The  high  sensitivity  of  viscoelastic  materials  causes 
the  matrix  values  of  Equation  (2.2)  to  change  with  only  mild 
variations  in  excitation  and  environmental  conditions.  This 
results  in  finite  element  solutions  that  are  only  applicable 
over  a  narrow  range  of  frequency  or  temperature.  A  solution 
to  Equations  (2.2)  and  (2.3)  which  reflects  varying  matrix 
values  may  thus  be  obtained  by  the  segmentation  of  the 
equations  into  incremental  modal  frequency  response  finite 
element  analyses  steps.  This  composite  approach  gives  good 
results  but  in  practical  applications  may  become 
prohibitively  expensive  and  time  consuming,  especially  when 
several  design  alternatives  are  under  consideration. 

An  application  of  the  modal  strain  energy  method  (MSE) 
by  Johnson  and  Kienholz  [Ref.  7s pp.  5-12]  is  suggested  as  an 
alternative  solution  to  Equations  (2.2)  and  (2.3)  by 
imposing  frequency  dependency  on  the  frequency  response 
through  post  analysis  data  manipulation.  In  the  MSE 
approach,  the  MSC/NASTRAN  code  is  used  to  construct  a  finite 
element  model  with  constant,  midrange  values  of  loss  factor 
and  shear  modulus  of  elasticity.  NASTRAN  modal  analysis  is 
applied  and  the  mode  shapes  and  natural  frequencies  are 
extracted  and  inserted  into  Equations  (2.2)  and  (2.3).   The 


fraction   of   elastic   strain   energy   attributed   to   the 

viscoelastic  core  at  each  natural  frequency  is  also  obtained 

from  the  modal   analysis   and  the  damping  material   loss 

factors  interpolated  at  each  natural  frequency.   The  matrix 

[2^{jOn]  of  Equation  (2.2)  is  then  replaced  with  [n(y— )(jOn] 

T 
when  Equations   (2.2)   and   (2.3)   are  modified  using  the 

interpolated  loss  factors,  the  frequency  response  solution 

becomes  frequency  dependent.    The  MSE  method  thus  gives 

accurate   solutions  over  a   large   frequency  range  while 

remaining  cost  effective. 
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Figure  2.1 


Variation  of  Shear  Modulus  and  Loss 
Factor  with  Temperature  [Ref.  5:p.  90] 


III.   MATERIAL  AND  METHODS 

A.  MODEL  DESCRIPTION 

Three  different  plate  structures  were  analyzed  during 
this  study  in  order  to  provide  a  basis  for  the  comparison 
with  the  finite  element  results.  These  plate  structures 
consisted  of  a  simple  plate,  a  milled  plate  with  an 
unrestricted  constraining  layer  and  finally  a  milled  plate 
with  a  restricted  constraining  layer.  The  simple  plate  was 
a  14  X  30  inch  aluminum  rectangular  plate,  3/8  inch  thick 
(Figure  3.1).  Dimensions  were  chosen  to  give  approximately 
15  resonant  frequencies  under  1000  Hz.  The  milled  base 
plate  and  its  dimensions  are  shown  in  Figure  3.2.  All 
plates  were  machined  from  standard  6061  T6  aluminum.  The 
use  of  a  damped  plate  structure  in  a  corrosive  environment 
requires  design  modifications  to  protect  the  viscoelastic 
material.  The  welded  plate  structure  of  Figure  3.2  is  such 
a  design  using  the  same  material  and  overall  thickness  as 
the  simple  plate.  This  milled  plate  was  subjected  to 
vibration  analysis  before  and  after  welding  to  detect 
reduced  damping  due  to  the  welding  process. 

B.  DESIGN  OF  VISCOELASTIC  AND  CONSTRAINING  LAYER 
THICKNESSES 

To  approximate  system  loss  factor,  the  plate  structure 

was  represented  by  a  three  layer  beam  system  and  subjected 


to  the  Bernoulli-Euler  type  analysis  described  by  M.L.  Drake 
and  developed  by  Ross,  Kerwin,  and  Ungar  [Ref.  8: pp.  1-6; 
Ref.  9].  The  viscoelastic  material  and  constraining  layers 
were  varied  incrementally  in  the  Bernoulli-Euler  equations 
(3.1)  to  (3.3)  below,  using  the  FORTRAN  program  listed  in 
Appendix  Ao 

3      3      3 


2 
+  E3H3 (H3^-D) 2  -  [:g-  +^(H2i-D)  +E3H3  (H3,-D)  ]'^  (3.1) 


E^H^   E2H2  H^,-D 


^1 

_  E2H2  (H2-L  -  -2-)   +   g^E2H2H21  "^  ^3"3"31^  ,^  ,, 

° Eli; ^^'^^ 

Vl  '*"^T^  +g(E;L^]_  +E2H2  +E^H  ) 


^1  "^^  ^1  "'■^^  "^^3             S 

H^T  =  9 — ;   H-.,  =  5 ;  and  g  =  5-         (3.3) 

-^      ^  -^-^       ^               E3H^H2K^ 
where : 


E  =  Young's  modulus  of  elasticity. 

G  =  shear  modulus. 

I  =  second  moment  of  area. 

H  =  thickness. 

~  r  ^Beam  j^proximation  Lengths  2      ,Beam  Approximation  Widths  2, 

K^   =  ^^  Plate  Length  '     ^  Plate  Width  ^    ^ 

and 
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Subscript  1  refers  to  the  base  structure. 
Subscript  2  refers  to  the  damping  layer. 
Subscript  3  refers  to  the  outer  layer. 
No  subscript  refers  to  the  composite  system. 

The  maximum  system  loss  factor  and  the  temperature  at  which 
it  occurred  were  found  for  each  step.  The  results  of  70 
incremental  thickness  changes  were  plotted  in  Figure  3.3. 

The  University  of  Dayton  Research  Institute,  in 
suggesting  damping  design  procedures,  recommends  a  system 
loss  factor  of  0.1  for  most  applications.  [Ref.  10] 
Because  Equations  (3.1) -(3. 3)  are  only  an  approximation,  the 
desired  system  loss  factor  was  increased  to  0.3  in  this 
analysis.  From  this  figure  and  based  on  a  loss  factor  of 
. 3 ,  the  values  of  the  viscoelastic  damping  layer  and 
constraining  layer  were  chosen  to  be  .0625  and  .25  inches 
respectively . 

C.   VISCOELASTIC  MATERIAL 

The  damping  material  applied  to  all  three  models  tested 
was  3M  SJ-201  ISD-112  viscoelastic  compound  with  a  thickness 
of  60  mil,  supplied  with  suitable  adhesive  already  applied 
to  both  sides  by  the  manufacturer,  3M,  St.  Paul,  Minnesota. 
During  application  of  the  viscoelastic  and  constraining 
layers,  care  was  taken  to  allow  air  to  escape  from  between 
the  layers.   Entrapped  air  can  result  in  severe  bond  defects 
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giving  lower  natural  frequencies  and  reduced  system  loss 
factor.   [Ref.  ll:pp.  8-9] 

ISD-112  was  chosen  for  its  excellent  loss  factor  and 
shear  modulus  characteristics  in  the  zero  to  100°  F 
temperature  range  for  frequencies  below  2,000  Hz.  This 
trait  is  exhibited  by  the  broad  loss  factor  curve  of  Figure 

3.4  and  is  highly  desirable  in  many  applications.   Figure 

3.5  shows  a  University  of  Dayton  nomogram  of  the 
experimental  determined  shear  modulus  and  loss  factor 
properties  for  ISD-112 .  To  determine  shear  modulus  and  loss 
factor  values,  locate  the  material  temperature  on  the  top  of 
Figure  3.5  and  the  frequency  of  interest  on  the  right 
margin.  The  diagonal  line  of  temperature  and  horizontal 
line  of  frequency  intersect  on  the  vertical  line  of  reduced 
frequency.  Reduced  frequency  is  the  frequency  value  of 
interest  multiplied  by  the  temperature  of  interest  and 
divided  by  the  reference  temperature  of  Figure  3.5.  Reduced 
frequency  and  reference  temperature  are  used  to  convert  the 
input  values  of  temperature  and  frequency  to  the  conditions 
used  to  plot  the  loss  factor  and  shear  modulus  curves. 
[Ref.  12:pp.  47-50]  The  intersection  of  reduced  frequency 
and  the  plotted  curves  give  the  horizontal  line  of  loss 
factor  and  shear  modulus  shown  in  the  left  margin.  While 
nomograms  are  useful  in  comparing  characteristics  when 
selecting  a  damping  material,  precise  values  are  difficult 
to  read  on  the  log-log  scale  and  the  actual  equations. of  the 
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nomogram  curves  were  used  to  give  accurate  values  of  loss 
factor  and  shear  modulus  in  the  finite  element  codes.  The 
equations  for  ISD-112  are  given  below  and  used  in  the 
program  listed  in  Appendix  A.   [Ref.  9: pp.  6-14] 

Log(FR)  =  log(F)  -  12* (T-TZERO)/ (525*T-TZERO)      (3.4) 

Log(G)  =  Log (ML)  +  (2*Log  10 (MROM/ML) )/ (1+FROM/FR) N    (3.5) 

Log(ETA)  =  Log(ETAFRO)  +  A*B+D* (1-SQRT (l+A^) )/ (C/2)     (3.6) 

where  the  ISD-112  viscoelastic  material  nomogram  constants 
are: 

TZERO  =131 

FROM  =  400,000 

MROM  =  2059.5 

N  =  .27 

ML  =  28.427 

ETAFRO  =  1.25 

SL  =  0.27 

SH  =  -0.3  5 

FROL  =  150,000 

C  =  1.3 
and: 

FR  =  reduced  frequency 

F  =  frequency 
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T  =  temperature 

G  =  shear  modulus 

ETA  =  loss  factor 

A  =  (Log(FR)  -  Log(FROL))/C 

B  =  SL+SH 

D  =  SL-SH 

D.   NASTRAN  FINITE  ELEMENT  MODEL 

Computer  analysis  of  the  sandwich  structure  was 
performed  using  the  MSC/NASTRAN  finite  element  program.  As 
shown  in  Figures  3.6-3.8  the  base  plate  and  constraining 
layer  were  constructed  of  QUAD4  plate  elements  with  3/16  and 
1/8  inch  offsets,  respectively.  The  more  complicated  HEXA8 
solid  elements  used  to  simulate  the  sandwiched  viscoelastic 
layer  were  required  to  accurately  represent  the  strain 
energy  due  to  shearing  of  the  viscoelastic  core  and  to 
account  for  the  coupling  between  stretching  and  bending 
deformations.  [Ref.  7:p.  8;  Ref.  13:Sec.  2.22,4.3]  The 
support  bulk  data  card  was  used  to  give  free- free  end 
conditions  to  the  plates.  After  constructing  the  finite 
element  models,  two  different  analyses  were  performed  to 
predict  frequency  response  amplitude. 

1.   Modal  Strain  Energy  (MSE)  Method 

Rearranging  Equations  (2.2)  and  (2.3)  and  solving 
for  the  response  amplitude  gives  the  following  equation: 
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^  {1  -  (oj/w  )^  ^  +  (n.  (^),  (co/oo^  J}^  ^ 

N    f^i<^>i<"/"n,i)><*E,i*R,i^*/"n^   -,  2\l/2 


-   I 


^  {1  -  (u/u  )^}^  +  (n.  (T^).(ui/ai„  h)}^ 
where :  -'- 


X(a3)  =  amplitude  of  response  as  a  function  of 
frequency 

^n  i  =  ith  natural  radian  frequency 

^E  i  ~  mode  shape  value  for  the  excitation  grid 
point  at  the  ith  natural  frequency 

*^R  i  ~  mode  shape  value  for  the  response  grid 
point  at  the  ith  natural  frequency 

Tlj^   =  iterated  viscoelastic  material  loss  factor 

(Vv/Vrp)  j^  =  fraction  of  elastic  strain  in  the  visco- 
elastic material  at  each  model 

L    =  magnitude  of  the  load  applied  at  the 
excitation  grid  point. 


The  MSE  method  utilizes  the  MSC/NASTRAN  method  to  find  the 
constant  values  on  the  right  side  of  Equation  (3.7).  The 
response  amplitude  is  then  calculated  external  to  the  finite 
element  code.  The  mode  shapes,  natural  frequencies  and 
fractions  of  elastic  strain  energy  in  Equation  (3.7)  were 
determined  using  a  single  MSC/NASTRAN  normal  mode  analysis 
for  each  plate.  The  damping  material  properties  input  to 
the  MSC/NASTRAN  normal  analysis  must  be  constant  and  were 
determined  using  the  nomogram  Equations  (3. 4) -(3. 6)  for  a 
temperature  of  60°  F  and  the  mid-range  frequency  of  500  Hz. 
The  output  of  the  MSC/NASTRAN  normal  analysis  gave  the 
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natural  frequency,  mode  shapes,  and  the  fractions  of  strain 
energy  in  the  damping  material  at  each  associated  natural 
frequency.  The  simple  plate  MSC/NASTRAN  normal  analysis 
program  is  listed  in  Appendix  B. 

The  loss  factor  of  Equation  (3,7)  is  sensitive  to 
the  natural  frequency  and  several  MSC/NASTRAN  mode  analyses 
must  be  applied  to  compensate  for  errors  in  the  derived 
natural  frequencies.  Johnson  and  Kienholz  describe  an 
iterative  method  to  approximate  this  loss  factor  frequency 
sensitivity.  [Ref.  7: pp.  27-31]  The  iterative  method  was 
applied  by  first  plotting  the  damping  materials  shear 
modulus  as  a  function  of  frequency  using  the  nomogram 
Equations  (3.4)-(3.6).  The  damping  material  shear  modulus 
of  each  plate  was  then  varied  through  several  MSC/NASTRAN 
modal  analyses  and  the  resulting  change  in  each  natural 
frequency  was  plotted  on  the  same  curve.  The  intersection 
of  the  curves  gives  an  approximation  of  the  natural 
frequency  and  shear  modulus  for  each  mode.  Figure  3.9 
illustrates  the  iterate  technique  for  two  of  the  simple 
plate  modes.  The  loss  factor  value  was  then  found  from  this 
approximate  shear  modulus  using  the  nomogram  Equations 
(3.4)-(3.6).  When  this  iterative  method  was  first  applied 
by  Johnson  and  Kienholz,  they  noted  the  results  differed 
from  those  observed  in  frequency  response  methods.  Modes 
below  the  mid-range  frequency  show  too  much  damping  and 
those  above  show  too  little.   [Ref.  7:p.  12].   Johnson  and 
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Keinholz  suggested  multiplying  the  loss  factor  obtained  in 
the  described  iterative  method  by  (Gj^/Gj-gf )  ^/^ ,  where  G-j^  is 
the  viscoelastic  material  shear  modulus  at  the  ith  natural 
frequency  and  Gj-gf  is  the  core  shear  modulus  at  the  mid- 
range  frequency.  Once  the  matrices  of  natural  frequency, 
mode  shape,  strain  energy  fraction  and  loss  factor  were 
calculated  for  each  plate,  Equation  (3.7)  was  solved  and  the 
frequency  response  amplitude  was  plotted  as  a  function  of 
frequency. 

2 .   Frequency  Response  Method 

Use  of  the  NMSC/NASTRAN  model  frequency  response 
analysis  to  directly  calculate  response  amplitude  within  the 
finite  element  program  requires  the  viscoelastic  material 
properties  to  be  input  as  constants  throughout  the  frequency 
range  analyzed.  This  restriction  limits  the  modal  frequency 
response  analysis  to  frequency  spans  over  which  the  visco- 
elastic materials  shear  modulus  and  loss  factor  have  little 
change.  Inspection  of  the  nomogram  in  Figure  3.5  indicated 
a  frequency  span  of  no  greater  than  100  Hz  would  give  little 
change  in  the  viscoelastic  material  properties.  In  order  to 
give  the  best  accuracy  at  the  peak  response  values,  the 
natural  frequencies  of  each  plate  were  estimated  by 
conducting  a  MSC/NASTRAN  modal  analysis  using  the  mid-range 
viscoelastic  material  properties.  The  viscoelastic  material 
shear  modulus  and  loss  factor  were  then  obtained  for  the 
plate  structures  at  each  of  these  natural  frequencies. 
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Modal  frequency  response  analysis  was  then  performed  over 
100  Hz  spans  centered  at  each  natural  frequency.  The  output 
magnitude  versus  frequency  curves  for  each  resonant 
frequency  were  compiled  and  plotted  as  composite  curves. 
The  modal  frequency  response  analysis  in  MSC/NASTRAN  for  the 
simple  plate  is  provided  in  Appendix  B» 

E.   VIBRATION  GENERATOR  TESTING  ARRANGEMENT 

Experimentally  determined  vibration  response  analysis 
utilized  a  Hllcoxon  research  F3  vibration  generator  as  the 
input  force  device.  The  F3  vibration  generator  was  attached 
to  the  base  plate  on  all  applications.  An  HP  3562A  Dynamic 
Signal  Analyzer  provided  the  swept  sine  driving  signal  which 
was  amplified  by  a  Wilcoxon  research  power  amplifier  (Model 
PA7C) .  A  force  transducer  was  located  before  the  attachment 
port  at  the  base  of  the  vibration  generator  to  measure  the 
input  excitation  force.  The  force  signal  was  amplified  by  a 
charge  amplifier  and  fed  as  the  input  signal  to  the  HP  3562A 
Dynamic  Signal  Analyzer.  An  Endevco  accelerometer  was 
placed  at  a  symmetrically  opposite  point  on  each  plate  and 
coupled  with  a  signal  conditioner  to  provide  the  output 
signal  to  the  HP  3562A  Dynamic  Signal  Analyzer.  Figures 
3.6-3.8  show  the  locations  of  the  vibration  generator  and 
accelerometer.  The  HP  3562A  Dynamic  Signal  Analyzer 
converted  the  signals  to  digitized  data  and  computed  the 
frequency  response,  coherence  and  power  spectra. 
Descriptive  points  were  taken  from  the  HP  3562A  Dynamic 
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signal  Analyzer  display  of  frequency  response  and  plotted  on 
the  same  graph  as  the  finite  element  predictions. 
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Figure  3.1   Dimensions  of  Simple  Plate  with  Applied 
Viscoelastic  and  Constraining  Layers 
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Figure  3.2   Dimensions  of  Milled  Plate  with  Applied 
Viscoelastic  and  Constraining  Layers 
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EIXCITATION  POINT 


1ZZP0HSZ    MEASURCMENT 


1/16"ISD-112  Viscoelastic 
Material  HEXA  Elements 


12   3   4. 


2,33  IN  Trpiau- 


3/8"  Aluminum  Base  Plate 
QUAD4  Elements 


Figure  3.6   MSC/NASTRAN  Finite  Element 
Model  for  Simple  Plate 
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EXCITATION  POINT 


1/4"  Aluminum 
Constraining  Layer 
QUA04  Elehients 


RESPONSE 
MEASUREMENT 


'•  ISO-112 
astic  Material 
:   F/ements 


Rim  of  Aluminum  Base  Plate 
(Upper  Portion)  Encircling 
Viscoelastic  Layer  and 
Constraining  Layer  3/8" 
QUA04  Elements 


1  IN  TYPICAL  - 

2  IN  TYPICAL- 


-**■ 


3/8"  Aluminum  Base  Plate 
QUAD4  Elements 


Figure  3,7   MSC/NASTRAN  Finite  Element  Model  for 

Milled  Plate  with  Free  Constraining  Layer 


26 


EXCITATION  POINT 


RESPONSE   MEASUREMENT 


2  IN  TYPICAL 


1/16"  ISD-112  Viscoelastic 
Material  HEXA  Elements 


4"  Aluminum 
training  Layer 
4*  Aluminum 
Base  Plate 
QUAD4  Elements 


1/16"  Aluminum  Base  Plate 
Encasing  Viscoelastic  Layer 
HEXA  Elements 


Aluminum  Base  Plate 
QUA04  Elements 
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Figure  3.8 


MSC/NASTRAN  Finite  Element  Model  for  Milled 
Plate  with  Welded  Constraining  Layer 
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IV.   RESULTS 

A  series  of  response  amplitude  versus  frequency  plots 
were  constructed  for  each  plate.  The  first  set  of  figures 
(Figures.  4.1-4.3)  plots  the  experimental  data  for  each 
plate  before  and  after  the  application  of  the  viscoelastic 
damping  and  constraining  layers.  After  completing  the  two 
prediction  methods,  the  results  were  plotted  (Figures  4.4- 
4.6)  with  the  experimental  analysis  for  the  damped  plates. 

A.   SIMPLE  PLATE 

The  experimental  frequency  response  plot  of  the  simple 
plate  (Figure  4.1)  shows  a  large  increase  in  damping  when 
the  viscoelastic  and  constraining  layers  are  applied.  All 
natural  frequencies  shift  lower  in  frequency  and  in  some 
cases  amplitude  peaks  are  combined  into  a  single  peak.  The 
system  loss  factor  for  the  first  mode  increased  from  2.2%  to 
34%  and  the  third  mode  increased  from  1.6%  to  14%  after  the 
application  of  the  viscoelastic  material  and  constraining 
layer. 

Comparison  of  the  two  simple  plate  prediction  methods 
with  each  other  and  the  experimental  results  show  several 
differences  (Figure  4.4).  The  MSE  curve  resonant 
frequencies  differ  from  the  NASTRAN  composite  curve,  with 
the  MSE  resonant  frequencies  being  slightly  higher  than  the 
NASTRAN  composite  curve  at  frequencies  below  500  Hz  and 
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slightly  lower  at  frequencies  above  500  Hz.  The  resonant 
frequencies  of  the  experimental  data  are  lower  than  either 
prediction  method  throughout  the  frequency  range  inspected. 
The  experimental  data  response  peak  amplitudes  are  higher 
than  either  of  the  prediction  curves.  The  response  peak 
amplitudes  of  the  MSE  method  have  approximately  the  same 
value  as  the  NASTRAN  composite  curve  with  the  poorest 
matching  above  600  Hz  and  the  best  matching  at  the  first 
peak.  The  experimental  data  show  a  loss  factor  of  34%  for 
the  first  mode  while  both  prediction  methods  indicate  a  45% 
loss  factor.  The  third  mode  has  a  14%  loss  factor  for  the 
experimental  data  and  a  20%  loss  factor  for  the  prediction 
methods . 

B.   MILLED  PLATE,  FREE  CONSTRAINING  LAYER 

The  experimental  frequency  response  plot  shows  a  large 
decrease  in  peak  response  amplitudes  after  the  viscoelastic 
damping  treatment  (Figure  4.2).  In  several  cases,  multiple 
amplitude  peaks  are  again  combined  into  a  single  peak.  The 
system  loss  factor  in  the  300  Hz  range  was  1%  prior  to 
damping.  Application  of  the  damping  material  increases  the 
system  loss  factor  to  4%  in  this  range. 

The  prediction  methods  were  plotted  for  the  finite 
element  model  in  Figure  4 . 5  with  the  two  prediction  methods 
showing  good  correlation.  At  frequencies  below  500  Hz  both 
prediction  methods  show  similar  peak  response  amplitudes  but 
at  differing  frequencies  with  the  NASTRAN  composite  method 
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having  consistently  lower  peak  frequencies.  Above  500  Hz 
both  prediction  methods  again  show  similar  amplitudes  but 
the  MSE  method  showed  lower  peak  frequencies.  The 
prediction  curves  are  poorly  matched  above  900  Hz.  The 
experimental  data  show  similar  resonant  frequencies  but  the 
frequency  peaks  are  sharper  than  either  prediction  method 
over  the  observed  frequency  range.  The  experimental  system 
loss  factor  at  the  3  00  Hz  peak  is  4%  while  both  prediction 
methods  indicate  a  6%  loss  factor. 

C.   MILLED  PLATE  WITH  WELDED  CONSTRAINING  LAYER 

After  the  experimental  analysis  of  the  unwelded  milled 
plate,  the  plate  constraining  layer  was  welded  to  the  base 
plate  while  maintaining  the  majority  of  the  plate  surfaces 
below  2  00°  F  during  the  welding  process.  After  welding,  the 
aluminum  constraining  layer  and  base  plate  were  visibly 
warped  and  the  plate  resonated  loudly  when  struck.  The 
experimental  analysis  of  the  welded  plate  confirmed  a 
significant  decrease  in  damping  (Figure  4.3).  The 
experimental  loss  factors  of  the  welded  plate  at  all  natural 
frequencies  were  indistinguishable  from  the  experimental 
loss  factors  for  the  milled  plate  prior  to  the  application 
of  the  viscoelastic  layer.  AFter  the  weld  was  cut,  the 
plate  once  again  exhibited  high  damping  (Figure  4.7).  The 
plot  of  the  two  finite  element  predictions  (Figure  4.6) 
yielded  curves  matched  for  frequencies  below  700  Hz  and 
curves  beginning  to  disagree  above  700  Hz. 
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V.   CONCLUSIONS 

Figure  4.1  shows  the  simple  plate  vibration  response  is 
reduced  significantly  after  the  3M  ISD-112  viscoelastic 
material  and  constraining  layer  are  applied.  Comparison  of 
the  experimental  data  and  the  two  prediction  methods  in 
Figure  4.4,  however,  show  that  both  the  NASTRAN  composite 
curve  and  the  MSE  method  predict  more  damping  than  has 
occurred.  Vibration  energy  absorbed  by  the  viscoelastic 
layer  is  converted  to  heat  and  possibly  increased  the 
damping  material  temperature.  No  attempt  to  hold  the  plate 
temperature  constant  during  the  experimental  analysis  was 
made.  Figure  3.9  shows  that  increasing  the  temperature  of 
the  viscoelastic  material  will  decrease  its  loss  factor. 
Another  contributing  factor  to  the  lower  actual  damping  may 
be  the  variation  in  the  behavior  of  the  viscoelastic 
material. 

In  examining  the  MSE  method  (Equation  (3.7)),  it  can  be 
seen  that  the  off-resonant  frequency  terms  have  an 
increasing  contribution  to  the  predicted  response  amplitude 
as  the  frequency  increases. 

The  natural  frequencies  of  the  MSE  method  are  found 
using  constant  mid-range  frequency  viscoelastic  material 
properties.  This  results  in  the  MSE  method  amplitude  peaks 
of  Figure  4.4  shifting  downward  from  the  NASTRAN  composite 
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curves  in  the  high  frequency  range.  An  upward  shift  of  the 
NASTRAN  curves  occurs  in  the  low  frequency  range.  The  two 
prediction  methods  for  the  simple  plate  show  more  disagree- 
ment than  for  the  other  two  plate  models  investigated.  The 
MSE  summation  Equation  (3.7)  contained  18  modes  for  each 
plate  analysis  conducted.  However  the  simple  plate  is  the 
least  stiff  and  thus  14  of  these  modes  occur  below  1000  Hz 
(78%) .  The  milled  plate  had  10  modes  below  1000  Hz  before 
welding  (56%)  and  8  after  welding  (44%)  .  The  number  of 
terms  in  the  simple  plate  summation  of  Equation  (3.8)  should 
contain  more  terms  to  increase  the  accuracy  of  the  MSE 
method  prediction  curve. 

Inspection  of  the  experimental  analysis  in  Figure  4.2 
again  shows  significant  damping  after  application  of  the 
viscoelastic  damping  material  and  constraining  layer. 
Figure  4 . 5  shows  that  both  predicted  damping  levels  were 
greater  than  the  experimental  data  with  the  causes  of  the 
error  being  the  same  as  for  the  simple  plate.  Examination 
of  the  two  prediction  curves  shows  the  MSE  method  gives 
acceptable  values  of  both  amplitude  and  peak  response 
frequency  throughout  the  frequency  range  analyzed.  The  MSE 
prediction  method  of  the  milled  plate  with  a  free 
constraining  layer  again  shows  the  MSE  method  response 
frequency  shifts  described  in  the  simple  plate  MSE 
prediction. 


40 


After  welding  the  constraining  layer  to  the  base  plate, 
the  structure  was  warped  and  exhibited  no  damping 
contributable  to  the  viscoelastic  material  layer  (Figure 
4.6).  When  the  weld  was  cut,  the  plate  warpage  decreased 
significantly  and  the  milled  plate  once  again  displayed  a 
large  loss  factor  (Figure  4.7).  This  indicated  that  the 
viscoelastic  layer  had  separated  from  the  plates  as  a  result 
of  the  warpage  caused  by  the  edge  welding  process.  Because 
of  the  loss  of  damping  through  deformation  of  the  plates 
during  welding,  only  the  two  prediction  methods  for  the 
milled  plate  with  a  welded  constraining  layer  are  plotted  in 
Figure  4.6.  This  plate  was  the  stiff est  plate  analyzed, 
having  eight  natural  frequencies  in  the  frequency  range 
studied  vice  10  natural  frequencies  for  the  milled  plate 
prior  to  welding  and  14  for  the  simple  plate.  As  a  result, 
the  response  amplitude  of  the  MSE  prediction  method  is  in 
agreement  with  the  NASTRAN  composite  prediction  method.  The 
MSE  prediction  method  again  shows  a  slight  upward  shift  of 
the  low  frequency  range  peaks  and  a  downward  shift  of  the 
high  range  peaks. 

Overall  indications  of  the  three  plates  show  actual 
plate  behavior  is  accurately  estimated  by  both  prediction 
methods.  Errors  occurring  in  the  higher  natural  frequencies 
are  due  to  an  insufficient  number  of  terms  in  the  MSE 
summation  Equation  (3.4).  Where  inaccuracies  occur,  the  MSE 
prediction  method  is  found  to  be  a  good  approximation  of  the 
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NASTRAN  prediction  method  and  should  be  an  acceptable 
substitute  for  the  NASTRAN  composite  prediction  method.  The 
simplicity  of  application  coupled  with  the  low  computer  time 
required  make  the  MSE  method  the  most  attractive  choice  for 
engineering  analysis  which  must  consider  both  a  wide  range 
of  environmental  conditions  and  several  design  configurations. 
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VI.   RECOMMENDATIONS 

The  following  are  suggestions  for  further  studies  in 
this  field: 

1.  This  study  investigated  prediction  accuracies  of  the 
MSE  and  composite  method  at  constant  temperature. 
Accuracy  over  a  temperature  span  should  be 
investigated . 

2 .  The  results  of  the  MSE  method  were  heavily  dependent 
on  the  MSC/NASTRAN  normal  mode  analysis  conducted 
with  viscoelastic  properties  evaluated  at  500  Hz.  An 
investigation  of  the  prediction  sensitivity  to  shifts 
in  evaluation  frequency  is  recommended. 

3.  The  use  of  the  HEXA8  solid  finite  element  to  model  the 
viscoelastic  layer  gives  results  in  an  MSC/NASTRAN 
program  which  requires  significant  computer  time  to 
solve.  Use  of  the  viscoelastic  beam  element  to  model 
the  viscoelastic  layer  have  been  suggested  to  reduce 
computer  time  requirements.   [Ref.  14: pp.  352-354] 

4 .  Failure  of  the  welded  plates  viscoelastic  material  was 
due  to  the  high  temperature  gradient  causing  warpage, 
yet  failure  to  cool  the  surfaces  of  the  plate  would 
have  resulted  in  destruction  of  the  damping  material. 
Investigation  of  welding  techniques  is  required  to 
give  practicality  to  the  milled  plate  design  used  in 
this  paper. 

5.  Application  of  several  viscoelastic  material  and 
constraining  layers  have  been  used  on  structures  in 
order  to  broaden  the  frequency  range  at  which  damping 
occurs.  Applying  the  prediction  methods  investigated 
to  such  multiple  damping  layers  should  be  researched. 

6.  Viscoelastic  material  properties  may  vary  from  batch 
to  batch  of  the  same  material.  The  construction  of  a 
nomogram  using  samples  of  the  batch  applied  to  the 
structure  should  greatly  increase  prediction 
accuracies. 
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APPENDIX  A 
NOMOGRAM  LOSS  FACTOR  AND  SHEAR  MODULUS  CALCULATIONS 

As  explained  in  Section  III. A,  it  is  impractical  to  use 
nomogram  chart  values  of  loss  factor  and  shear  modulus  in 
analyzing  damping  characteristics  of  systems  containing 
viscoelastic  materials.  All  calculations  of  shear  modulus 
and  loss  factor  in  this  paper  were  accomplished  through  use 
of  the  FORTRAN  program  listed  in  Figures  A.1-A.4.  Although 
this  program  is  not  complicated,  use  of  the  equations 
without  computer  assistance  is  extremely  tedious  due  to  the 
large  number  of  terms  e  Note  that  page  4  of  the  program 
(Figure  A. 4)  also  gives  the  system  loss  factor  based  on  the 
Bernoulli-Euler  equations  (3.2)-(3.3).  Equations  (3.2)  and 
(3.3)  require  Young ' s  Modulus  and  the  shear  modulus  values 
in  order  to  be  solved  and  thus  it  is  most  convenient  to 
include  these  calculations  here. 
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APPENDIX  B 

MSC/NASTRAN  CODE  FOR  SIMPLE  PLATE  NORMAL  MODE  ANALYSIS 
AND  MSC/NASTRAN  CODE  FOR  SIMPLE  PLATE  FREQUENCY 

RESPONSE  ANALYSIS 


While  the  MSC/NASTRAN  code  was  used  for  all  the  finite 
element  plates  analyzed,  listing  each  program  would  serve  no 
useful  purpose.  For  completeness,  a  normal  mode  analysis 
and  a  frequency  response  analysis  are  listed  here  for  the 
simple  plate  structure  (Figures  B.l-B.lO).  The-  MSC/NASTRAN 
programs  for  the  milled  plate  with  free  constraining  layer 
and  milled  plate  with  welded  constraining  layer  vary  only  in 
the  geometric  organization  of  the  grid  points  and  finite 
elements. 
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ID  DY^<AMICS#FREQ  MOD 

SDL  3 

TIME  2^0 

DIAG  8f9,13 

CEND 

TITLEsMODAL  A^JALYSIS  500  HZ  PROPERTIES 

ECHOsUNSORT 

SPCslo 

METHODsl 

ESE=ALL 

DISP=ALL 

OJTPUT(PLOT) 

PLOTTER  NASTRAN 

PAPER  SIZE  U.  BY  10. 

A)(ES  =  Z,XrY 

SET  1=ALL 

FIND, SCALE, SET  1, ORIGIN  1 

PLOT 

PLOT  MODAL'  DEFORMATION  0, 1 , THRU, 15, SET  1, ORIGIN  1, SHAPE 

BEGIN  BULK 

GRID, 1,0,0,0,0 

2,*(7),s,s,*(S.75),s 

=  7 

GRID, 2, 0,2. 5333, 0,0 

=,*(7),=,s,*(3.75),= 

=  7 

GRID, 5,0, a, 6667, 0,0 

=,*(7),=,=,*(3.75),2 

=  7 

GRID, a, 0,7. ,0,0 

=,*(7),=,=,*(5.75),s 

s7 

GRID, 5^0, 9. 5333*0,0 

s,*(7),2,s,*(3«75),s 

s7 
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10  DYMAMlCSfFREQ  MOD 

SDL  30 

TIME  2ao 

CEND 

TITLE=FREauENCY  RESPONSE  500  HZ  PROPERTIES 

ECHO=JNSORT 

SPC=10 

METH0D=1 

Sy/ECTOR  =  ALL 

OLOADs2000 

FREQUENCY=2100 

SET  lll=2a 

DISPLACEMENT (PHASE) =111 

BE5IN  BJLK 

GRID, 1,0,0,0,0 

=,*(7),r,=,*(3.75),s 

=  7 

GRID,2,0,2,3333,0,0 

s,*(7),s,r,*(3.75),= 

=  7 

GRID, 3#0, a. 6667, 0,0 

s,*(7),=,=,*(3.75),s 

=  7 

GRID, a, 0,7. ,0,0 

r,*(7),=,=,*(5,75),x 

=  7 

GRID, 5, 0,9, 3333, 0,0 

=,*(7),=,=,*(3.75),= 

=  7 

GRID, 6,0, 1 1.6667,0,0 

=,*(7),r,=,*(3.75),= 
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GRID,7,0,ia.,0,0 

r,*(7),=,=,*(3.75),= 

=  7 

GRID,  101,0,0,0, .0625 

s,*(7),s,=,*(3.75),s 

=  7 

GRID, 102, 0,2. 3333,0,. 0625 

=#*(7),=,s,*(3.75),s 

s7 

GRID, 10 3,0, a. 6667,0,. 0625 

=,ft(7),=,s,*(3.75),= 

=  7 
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GRID,lOa,0,7,,0, .0625 

s,*(7),=r=,*(3.75),= 

=  7 
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=,*(7),=r=,*(5,75),s 
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CaUADa,l,10  0,l,2,9,8,,-.1875 


=,*(6),=,*(7),*(7),*(7 

=  6 

CaUADa,2, 100,2,3, 10,9, 

=,*(6),=,*(7),*(7),*(7 
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=,*(6),=,*(7),*(7),*(7 
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C3UA0a,6, 100,6,7, la, 13 

=,*(6),=,*(7),*(7),*(7 
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3,»(6),=,»(7),*(7),*(7 
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CQUA04,102,101,102,103 

=,*(6),=,»(7),*(7),*(7 

=  6 

C3JADa,103,101,103,10a 
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,-.1875 
#*(7),,= 
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